«TvAL  POSTGKADUAT*  SO* 


sicticw 

JATE 
WOHTKHET.CAIJFOaNiA 


NAVAL 


*»*C 


NP£53=7B-001 


NAVAL  POSTGRADUATE  SCHOOL 


? 


Monterey,  California 


REPRODUCING  KERNEL  FUNCTIONS  FOR  THE 
SARD  CORNER  SPACES  B 


fp>ql 


RICHARD  FRANKE 


Interim  Report  for  Period  January  -  March  1978 


Approved  for  public  release;  distribution  unlimited 


ired  for:  Chief  of  Naval  Research 
feddocs  Arlington,  VA  22217 

D  208.14/2:NPS-53-78-001 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  Tyler  F.  Dedman  Jack  R.  Borsting 

Superintendent  Provost 


The  work  reported  herein  was  supported  by  the  Foundation  Research 
Program  of  the  Naval  Postgraduate  School  with  funds  provided  by  the 
Chief  of  Naval  Research. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 

This  report  was  prepared  by: 


UNCLASSIFIED 


SECURITY   CLASSIFICATION   OF   THIS   PAGE  (When   Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.     REPORT   NUMBER 

NPS53-78-001 


2.  GOVT    ACCESSION   NO 


3.     RECIPIENT'S  CATALOG   NUMBER 


4.     TITLE  (and  Subtitle) 


Reproducing  Kernel  Functions  for  the 
Sard  Corner  Spaces  Br-  -j 


5.     TYPE  OF    REPORT  &   PERIOD  COVERED 


Interim  January  -  March  1978 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!"*.) 


8.  CONTRACT  OR  GRANT  NUMBERS 


Richard   Franke 


9.     PERFORMING  ORGANIZATION   NAMF.   AND   ADDRESS 


Naval  Postgraduate  School 
Monterey,  CA  93940 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a   WORK  UNIT  NUMBERS 

61152N  ,  RR 000-91-0]  , 
N0001478WR80023 


II.     CONTROLLING  OFFICE  NAME   AND  ADDRESS 

Chief  of  Naval  Research 
Arlington,  VA  22217 


12.     REPORT  DATE 

March   1978 


O.     NUMBER  OF  PAGES 

8 


14.     MONITORING  AGENCY  NAME  4   AODRESSC//  dlfterent  from  Controlling  Office) 


13.     SECURITY  CLASS,  (of  thta  report) 


UNCLASSIFIED 


15*.     DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (of  the  eb, tract  entered  In  Block  20,   It  different  from  Report) 


IB.     SUPPLEMENTARY   NOTES 


19.     KEY   WORDS  (Continue  on  reverae   elde  If  necaaaary   and  Identify  by  block  number) 


approximation  theory,  Sard  corner  spaces,  reproducing  kernel  functions 


20.     ABSTRACT  (Continue  on  reverae  aide  It  neceaaary  and  Identity  by  block  number) 

The  reproducing  kernel  functions  for  the  Sard  corner  spaces  B     are 

Tp>ql 

integrated  to  obtain  explicit  expressions.  The  smoothness  properties  of 
these  functions  are  discussed. 


DD     1   j*N   73     1473  EDITION  OF    1  NOV  65  IS  OBSOLETE 

S/N    0  102-014- 6601   | 


UNCLASSIFIED 


SECURITY  CLASS'FICAT'OM  OF  -uss  PAGE  (Whan   Data  Entered) 


1 .   Introduction 

Reproducing  kernel  functions  arise  naturally  in  approximation  theory  in 
Hilbert  spaces,  and  are  a  useful  theoretical  tool,  and  can  also  be  used  for 
computations.  Since  I  have  used  reproducing  kernels  in  Sard  corner  spaces 
in  my  own  work,  I  have  felt  it  would  be  useful  to  have  the  formulas  (given 
in  terms  of  integrals)  evaluated  and  written  out  in  a  more  explicit  form. 
Therefore  this  small  report  was  written  to  discuss  Sard  corner  spaces  and  to 
develop  those  formulas. 


2.  Sard  Spaces 

Sard  spaces  of  type  B  are  spaces  of  functions  which  have  a  certain 
type  of  Taylor  series  expansion  with  remainder  [2].  The  partial  derivatives 
which  exist,  and  in  terms  of  which  the  expansion  is  given,  is  specified  by 
the  "complete  core".  For  details  the  reference  should  be  consulted.  We 
follow  the  notation  of  Sard  so  that  subscripts  denote  differentiation.  The 
complete  core  of  Br  -i  ("  B  corner  p,q  ")  is  made  up  of  the  elements 

"S,t=   {fp,q(5't)} 


3s,b 


;a,t 


=  {fn  ,(s,b):  j  <  q} 
=  {f.  n(a,s):  i  <  p} 


w_  .  =  {f.  ,(a,b):  i  <  p,  j  <  q} 


Here  (a,b)  is  a  fixed  point. 

The  differentiation  diagram  is  at  right. 


argument  (a,t) 


argument  (a,b) 


*    t 


(p.q) 


(0,0) 


argument  (s,b) 


-1 


A  function  fe  Br  -i  can  be  expressed  as 

f(s.t)  =  I   (s  -  a)(i)(t  -  b)(j)f.  .(a,b) 
i<p  1,J 

j<q 

*    lis-   a)(i)/  (t  -  t)(q-1}f.  fl(a,t)dt 
i<p        b  1,q 

+  I(t-  b)(j)/S(s  -?)(P-])f  .(s,b)ds 
j<q        a  P'J 

+  /S/t(s  -s^Ct-  t)(q-1}f  fl(s,t)dtds 

a  b  p,q 

The  inner  product  for  Br  -i  ,  as  defined  by  Barnhill  and  Nielson  [1],  is 


(f,g)  =  I   f..(a,b)g.  (a,b) 
i<p  1J     1,J 


■P 

j<q 

6 
+  I    I  U   a(a»t)g.  (a,t)dt 
i<p  6  1,q     1,q 

+  I  Ap  j(s.b)g  ,(s,b)ds 

j<q  a  P'J      P'J 

/°7  f  q(s,t)g   (s.t)dtds 

a  6  P'q      P,q 

Here  [a, a]  x  [g,8]  is  the  rectangle  which  contains  the  region  of  interest 

The  above  leads  to  the  reproducing  kernel  function,  which  in  turn  will 

allow  construction  of  the  representers  of  functionals,  particularly  point 

evaluation  functionals,  in  which  we  are  interested.  The  reproducing  kernel 

function  is 

(*)  K(u,v;s,t)  =  I  (u   -  a)(1)(v  -  b)(j)(s  -  a)(i)(t  -  b)(j) 
i<p 

j<q 

+  I  (u  -  a)(i)(s  -  a)(i)/V(v  -  t)(q-])(t  -  t)(q-1)!|,(b,t,t)dt 
i<p  b 
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+    J  (v  -  b)(j)(t  -  b)(j)/U(u  -  s)(P-^(s  -  s)(p-1)^(a,s,s)di' 

j<q  a 


u  v 


+  /  /  (u  -  s)(P-])(s  -  s)^-1^  -  tJ^^Ct  -  t)(<H^(b,t,tMa,s,s)dtd? 
a  b 

In  the  above  the  function     iji     is  defined  as 

1       if     a  <  s  <  s 
4>(a,s,s)  {    -1       if     s  <  s   <  a 

0       otherwise 

For  the  above  case,  the  kernel   function  can  be  factored  into 
K(u,v;s,t)  =  g  (a;u,s)  g  (b;v,t),  where 

g  (a;u,s)     =     I  (u  -  a)(i  ]  (s  -  a)(i }  +  /"(u  -  sVP~]  }(s  -  s) (p_1  ^(a,s  ,s)ds 
p  i<p  a 

This  fact  simplifies  investigation  of  the  properties  of  the  representers  of 

functionals  considerably.  It  is  desirable  to  obtain  an  expression  for 

g  (a;u,s)  which  does  not  involve  the  integral.  Repeated  integration  by 

parts,  and  assuming  that  a  <  u,s  yields 

gp(a;u,s)  =  (-l)p(u  -  s)fp-]) 

+  I   {(u  -  a)(i)(s  -  a)(i)  +  (-l)1(s  -  a)^-1'1^  -  a)(p+i)}  . 
i<p 

Note  here  that  we  have  lost  the  obvious  symmetry  of  g   in  u  and  s  , 

although  not  the  actual  symmetry,  of  course.  We  can  regain  it  however,  at 

the  expense  of  two  formulas  for  g  ,  depending  on  whether  u  <  s  or  s  <  u 

To  do  this  expand  (u  -  s)|  p~  '  (for  u  >  s)  in  a  binomial  way,  obtaining 

f-nPfu       s^P"1)  -   ,H)P     2?"VD  n,/   n2p-l-k   (u-a)k  (s-a)2^1"* 
t-U    CU  -  sj  "   (2p-l)!    k^Q   (2p-l)!(-l)  k!      (IFT^T)! 

=   (-l)P2f1(-1)2P-1-k(u-a)(k)(s.a)(2p-l-k) 
k=0 


-3- 


=  (.1)Ppj1H)2p-i-k(u.a)(k)(s.,)(2p-l-W 

k=0 

M-D'T1!-!)2^-")1^*-*)^1 
k=P 


i<p 


I   (U  _  a)(p+i)(s  -  a)(p-1"1,M)1    (k*-p+i) 


i<p 


Thus 


,  (-l)p(u  -  s)(2P-])  +  I   (u  -  a)(P+i)(s  -  a)(P-1-'f)(-l)i 


i<p 


i<p 
and  we  see  that  if  we  interchange  u  and  s  the  expression  for  g   is 
unchanged,  and  thus  we  can  take 

(1)  gD(a;u,s)  =  I   {(u  -  a)(i)(s  -  a)(i)  +  (-l)1(u  -  a)(P+i)(s  -  a)^1"1*} 

H        i<p 

when  a  <  u  <  s  ,  and 

(2)  g  (a;u,s)  =  g  (a;s,u)  when 

a  <  s  <  u  .  When  s  <  a  <  u  we  find  that 

(3)  gD(a;u,s)  =  [  (u  -  a)(i)(s  -  a)(i) 

p  i<p 

For  u  <  a  it  is  easily  observed  that  the  only  difference  is  that  the 

function  ^(a,s,s)  changes  sign,  hence  the  terms  from  the  integral  appear 

with  opposite  sign,  so  we  can  say 

(T)  g  (a;u,s)  =  [  {(u  -  a)(i)(s  -  a)(i)  -  (-1)%  -  a)(p+1)(s  -  a)(p-M)> 
p        i<p 

when  s  <  u  <  a  ,  and 

(21)  g  (a;u,s)  =  g  (a;s,u)  when  u  <  s  <  a  .  When  u  <  a  <  s  we  find  that 

r  r 

(3')  g_(a;u,s)  =  £  (u  -  a)(i)(s  -  a)(i)  . 
H        i<p 
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It  is  possible  to  write  the  reproducing  kernel  function  in  a  somewhat 
more  concise  manner  by  defining  the  function 

Gp(a;u,s)  =  (-l)P(s  -  u)j2p_1) 

+  I{(u-a)(i)(s-a)(i)+(-l)i(u-a)(P-1-i)(s-a)|P+i)}  , 
i<p 

for  a  ^u  ,  and  then  observing  that  for  cases  1,  2,  and  3 

gp(a;u,s)  =  Gp(a;u,s)  , 


while  for  cases  T,  2',  and  3'  , 

gp(a;u,s)  = 
For  p  =  1  and  a  fixed  value  of  u  ,  graphs  of  g^u.s)  are  given  below, 


gn(a;u5s)  =  G  (-a;-u,-s) 


(u,l  +  u  -a) 


s  =  a  s  =  u 


(u,l  +  a  -  u) 


N 

(a,l) 

i 

i 

s  =  u 

s  =  a 

-5- 


For     p  =  2   ,  the  graph  of  the  function     g2(a;u,s)     is  given  below  for 
a  <  u.     For     a  >  u   ,  the  graph  is  flipped  about     s  =  a   ,  as   in  the  p  =  1   case, 


g1   =  1  +  (u  -  a)(s  -  a) 

+  (s  -  a)(u  -  a)(2)-(u  -  a 


(u.l  +  (u  -  a)(2)  +  l<u  -  a)3) 


1   +  (u  -  a)(s  -  a)  +  (u  -  a)(s  -  a) 
-  (s  -  a)<3> 


g     =  1   +  (u  -  a)(s  -  a 


'2 

We  note  that  g2  is  C2  in  s  except  at  s  =  a  ,  where  it  has  a 
continuous  first  derivative,  but  not  a  continuous  second  derivative.  In 


general  g   is  C 


2p"2  in  s  ,  except  at  s  =  a  ,  where  it  is  C 


p-1 
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